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This paper, an expansion and discussion of a series of lectures given 
by Er. F. K . ochcubel of the Technical Institute, Earnetadt, Gejraany, was 
undertaken by the writer as a thesis project while working for the degree 
of Master of Science in Aeronautical Engineering at the University of Mich- 
igan, Arji Arbor, Michigan. These lectures were given at the University of 
Michigan during the fall of 1951 while Dr. Scheubel was visiting the Uni- 
versity at the invitation of Dr. F. W. Eonlon of the Aeronautical Engineer- 
ing < Depart nent . 

Dr. Scheubel presented the material in six lectvires, the first three 
devoted to quasi-static stabilitv concepts, the Last three to their use in 
the solution of the equations of motion. Hie limited time prevented a de- 
tailed accounting of the assumptions involved as well as minute explanations 
of the approximations made. It is the purpose here to verify and expand ,on 
his presentation to a degree consistent with the limited scope of such a 
paper. 

Briefly, his material covered the following. His lectvires were re- 
stricted to symmetric or longitudinal motion. He developed the longitudi- 
nal equations of motion and the quasi-static stability criteria at equili- 
brium and at constant speed. He solved the equations using these quasi- 
static stability concepts for both the phugoid and short period modes for 
the stick-fixed case. An amplitude-phase relation for the two variables 
applicable to the modes was discussed. Ho then introduced the degree of 
freedom about the elevator hinge line and solved for the stick-free case. 
Finally the effects of an elevator impulse were discussed. 

This approach to dynamic longituiinal mction differs fi'om standard 
methods mainly in the quasl-stetlc stability concepts as regards their in- 
sertion in the sclutions tc the ecuatlons of motion. In this respect, the 
method may, as illustrated, be safely applied only to conventional air- 
craft, l.e., a rigid body, where the number of degrees of freedom and thus 
the complexity of solution is restricted. The analysis of the motion as- 
suming a phygold and a short period oscillation serves only to Illustrate 
the handling of the quasi-static concepts. The limitation of the method 
in respect to acceptable separation of the motion into these modes was 



broxi^'ht out by Dr. Scheubei and is discussed. The form of the solution of 
the complete set foregoing this separation is Included in this paper. 

The equations of motion as established by the standard approach of 
reference 1 are relega.ted to the flight path direction or wind axes In 
Older that the further discussion will bo based on common assumptions 
and approximations. Dr. Scheubei 's equations as developed for a system 
relative to the flight path v;lll then be compared terra by term with the 
standard equations. From this point onward, i.e., bsglrjiing with non- 
dimensionalizing the equations, the presentation is as Dr. Scheubei de- 
veloped it. 



KOtoKCIATUBE 



m = W/g 


mass of the aircraft 


Up 


initial velocities along x and l body 
axis 


u' , w 


perturbation velocities along x and z 
body axis 


q 


angular velocity about y axis (body or 
wind) 


X, Z, M 


air forces and moments 


e, t , C/', or Ae, A r ,A 


perturbation of attitude angle, flight 
path angle, and angle of attack 


^>1 

y 


initial attitude angle 

aircraft radius of gyration about y 
axis 


I = 

yy y 

V = + u? 

oil 


moment of inertia about y axis 
Initial flight path velocity 


AY 


velocity perturbation along flight path 


V = V + AY 
o 


velocity along flight path following per 
turbation 


L, D, M 


lift, drag and air moment 


W 


aircraft weight 


T , E 
o' 0 


Initial thrust and drag 
elevator deflection 


"l 

f 


lift coefficient for aircraft 
air density 


s 


wing area 


s 


aircraft drag coefficient 


u = AY/Y 


non-dimensional velocity perturbation, 
ratio of perturbation to total flight 
path velocity 
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t = 

e g ^ VS 


adopted standard time unit 


- g sV2 


nvass density (See Appendix A) 


6 = t/t 


non-dimensional time 


M 

0 


initial moment about y axis 




tail length, distance from aircraft 
c.g. to center of pressure of horizon- 
tal tail 


z :* R± il 


conjugate complex root of characteristic 
equation. B is real part, I is imaginary 
part 


= S/l^ 

y y 




ap &2 or ApAg 


constant coefficients of equations 


X 


position of aircraft center of gravity 
from most forward part of aircraft in 
percent of wing chords 


X 

n 


position of neutral point of aircraft 
from most forward part of aircraft in 
percent of wind chords 




angle of attack of horizontal tail 


= Lh 


lift coefficient of horizontal tail 


"h 


area of horizontal tall 


Co 


angular frequency of mode of motion 


T , T 
p’ r 


periods of phygoid (p) and rotary (r) 
modes of motion 


L 

P 


wave length of phygoid 


CD 


distance between the hinge line and cen- 
ter of gravity of the elevator 


m 

e 


mass of the elevator 

radius of gyration of the elevator 




hinge moment coefficient of elevator 


"h 


mean aerodynamic chord of the horizon- 
tal tall 



elevator ImpulBe 



H constant altitude reference 

o 

h altitude change 

Subscript ® Indicates steady state. values in the response discussion. 
Any other symbols used are either obvious or are locally defined for ease 
in handling equations. 
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II. Development of the Method. 

A. The equations of longitudinal motion. 

The standard development in a perfectly general way as accomplish- 
ed in reference 1 result in the equations of longitudinal motion of the form, 



A.l 


m(u' + Wj^q) = X^.u' 


' + Xv+Xq-mg 
w q 


cos 


A. 2 


m(w - Uj^q) = 


+ Zw+Zq-rag 


sin 


A. 3 


mi^q=I q=M 


u' + M w + M q 






y yy u 


w q 





These equations are relative to body axes and contain the follow- 
ing assumptions and approximations: 

(a) Initial symmetric steady motion is assumed. 

(b) The air reactions do not depend on the rates of change of the 
variables, U, V and W, or their integrals. 

(c) Second order and higher terms of the air reactions are neglected, 
i.e., only infinitesimal disturbances from an initial steady motion are 
treated. 

(d) The aircraft has a plane of symmetey and the steady motion about 
which the disturbances occur is symmetrical with regard to that plane. 

(e) The disturbance initially imposed on the system is vinenforced and 
the controls are locked. This rules out an initial couple, M^, and the vari- 
ation of M with S . 

The variables are seen to be u, w and q. These are Illustrated in Figiire 

1. The relation between the body sixes and the wind axes is shown in Figure 

2. The body axes are primed. Angles are measured positive counterclockwise 

from the horizontal reference for y and from the wind line for iX . 

In regard to these figures it is seen that, 

(a) Since q is a velocity and 6 a displacement, 6 = /qdtsAfX + Ai*” 



' 1 ‘- • 
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(b) is the steady state velocity and, ~ ^ 



Then there is the definition, V = V + AV 

' o 



also. 



^ d- (AV) 



dt 



dt 



(c) In essence u is now AV, and w does not exist. 

The inertia terms of the equations of motion parallel and perpendicular 
to the wind direction end about the center of gravity become, neglecting the 
inertia force due to linear acceleration perpendicular to the wind direction, 

d(AV) 



for A.l: 
for A. 2: 
for A.3: 



m 



m V 



77 



dt 
dt 



d ( - 'S' ) 

dt~ 



The weight components relative to the wind axes become, 
along x: mg sin 

along z: -mg cos 

and these are only affected by A . 

There remains only the air reaction derivatives. The force created 
on the tail due to q is the largest resulting from this disturbance s.nd 
from experience it is known to be small along the wind axes so X and Z 

q q 

are neglected. 
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In the standard method the left and drag contributed to the air reac- 
tion forces both along the chord and nornial to it. Also velocity pertur- 
bations existed along both axes. Using v:lnd axes, cnly the velocity per- 
turbation, AV, nov exists. Any sinking velocity along the z wind axis is 
merely a change in angle of attack. The air reactions along x consist 
only of the drag in its entirety, neglecting variations of thrust. The 

di^g varies with both and V. The air reactions along z consist only' 

' ) 
of the lift in its entirety and is affected by the same quantities. A 

pure rotary perturbation about the center of gravity changes bcth o< and 
^’'but only CX affects the linear forces. The same applies to the air mo- 
ments. The variations with ^ are found by wind tunnel tests and any rota- 
tion of the model about the center of gravity is considered pure cX . This 
does not hold for the case of the effect on M of both A cx and d”'" . These 
tvo perturbations make up 0, which, when multiplied by the tail length, 
gives the effective sinking velocity of the horizontal tail. This sinking 
velocity in turn is felt as a change in angle of attack of the tail. 

Thus the exterr.al forces for A.l become, 

(mg sin ) AT 



and equation A.l is, 



n 



d(^V) SD 



Sr . 



dt 



^ ^ AY + mg cos ^ A 



Fquation A. 2 is, 



m V A 0^ + AV + mg sin S' A S' 



equation A. 3 becomes, 






dt^ 






60^ “ §v — 6'Si 
The use of the wind axes then has the advantage of the aircraft's for- 
ward motion being along the x axis so that the lift and drag forces always 
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lia along these axes and resolution into normal and chord forces is not 
necessary. Fvirthernore, no resolution of the velocity is required to he 
made along the axis of perpendicularity vhich is the case if the forces 
are resolved along the body axes. The disadvantage is that these axes 
change continuously so that the raocents and products of inertia change. 
For small disturbe-nces from initial horizontal flight these discrepancies 

I 

are considered negligible. Also the thrust does not act along the wind 
axes for all aircraft and, for an individual aircraft, at all tines. 

The continuity of this paper and the advantage of physical "feel" 
that is had from the development of the equations of motion from initial 
suppoeiti'cn of vind axes can be best maintained by including such a de- 
velopment in its entirety at this point. This is the method used by Dr. 
Scheubel and appears to be in general use in Germany. See Eeference 2. 



\ 




With the aircraft in an attitude and vith the external forces as 
shovn in Figure 3, the external forces are equated to the time rate of change 
of momentum in the direction of the wind axis, i.o., along V, and one equa- 
tion of equilibrium follows. 

~ = T - D - W sin 3”' 

S dt 



1C 



The external forces are caused to chan'^e hY email changes or pertur- 
bations in the values of , V, 'Jf' end T . The forces in equilibrium 
then, following a small disturbance, are the initial values plus the in- 

4 

crements due to the changes. 



- ^ = T - r - W sin Jf- 
g dt o o o 



St 

5v 



Sd 

§v 



Sd 



jN - ~ ^ ~ W cos A ^ \ 






This is the result of a Taylor Series expansion of the quantities T, 5 
and W as functions of , V, ^ and ^ in vliich second order terms and a- 

bove arc neglected under the assumption that the changes are smell. Assuming 
an initial condition of equilibrium, it is seen that this equation is compar- 
able to equation A.l with one exception. Thus all the assumptions prevlouelv 
enumerated hold. The exception is the thrust term. It is easily seen that 
thrust is not affected by changes in 0< , ^ , or ^ . The stick-fixed case 
is considered here also. Thus all terms involving ^ are neglected. 

This first equation is now manipulated into a form necessaiy for sub- , 
Sequent solution. Initial equilibrium and assuming S' small gives. 



T - B - W sin S' = 0 
o o o 



and, cos S' = 1 



or 



V = L 



also, L = C, t P v^S and, B = i P V^S 

^ ’ T ^ ^ B 

Multiplying terms in T by and ^ , and terms in B by and ^ is 

o c 

V St ^ V SB SB 



W dV = T ^ S- 

g dt o T^ SV 



V ■ ^o B 5v 



_V ST ^ S In T . _V ^B _ S In B 
T.. 5v ~ S In V " B S? " S In V 



V cc^ 
and letting 



o( - L 

.aV 



= u 



there is, 

W ^ _ 

g dt ■ 



V d(^) 



S 



dt 



dt 



B 



f r"<j 



T 



o S in T S In B 



‘B nrv ■ nr?) 

Using the standard notation, r—r- = and dividing bv r — S, there is 

- o( id 



S( 



g .a)P- 



O 0^ 
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V 

g 



2 dUv) 



f> V^S 



.' 1 CX 



, In D 

"L In V 



■*■0 3 In 

D h In V ^ “ 
o 



Thie eqxiatior. is non-direneione ll*.ed in the ueual manner. The right 
eide of the equation is non-dimensional as it stands. Considering the left 
side, piulti inlying ninerator hy V/V and the denominator bv 



W 

3 p v^S 



2 V_ d^'/v) 

t d(- 



s 






W 

A tine unit, t , is adopted bv letting, - ^ 

^ B 



■5— — = 1 then, 



"s = 






:ws 






]/2 



gyv, 

2W 



3/2=/^ V- = — vhere the mass number , is given by, 



gpS^ gp - 

t =6 

s 

60 that, d(^) = d 'Z’ and the variable ilV is transformed as, 

g 

The first equation of notion f^en becomes. 






(’) 



• 



Thie mass number, or density factor, is 



= u. 



A.I. H-u-.C. 



i(X - ^ b* - C, (r 

L 0 0 



5 1 



in o 



In V 



o S in T\ 

- 5 nj-v> “ •• 



- -O 

Again referring to Figure 3 > the external forces perpendicular to the 

vind line are equated to the centrifugal force or mas^ times the centrifugal 

. W „ dO W -- a _ ,, V*. , W d « ^ , 

scceleration. -V— = - V— = L-W cos y since, - 1 tt = C, because 
g at g Gt ' 3 dt ' 

V ^ results in no chonge in velocity perpendiculer to the flight path and the 
linear acceleration in this direction is neglected. If a small disturbance 
is introduced, and, since lift does not vary with ^ and the etick-fiied case 
is assumed, 

- V L - W cos if ■ ‘ - W sin ^ .1 r . 

g dt o o ocA dV 

Vith the sane assumption as before, that of initial equilibrium, this 

equation is seen to be the seme as equation A . 2 and so the same assumptions 

and approxlraa.ticr.c are implied. Expanding the derivatives. 
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end 



sl<'Lf r^> ■ ^?7 



•L^f-I'f r=^'o^faRr-7r> 

S In C, 



1 ^ O 

"'L V °'^L ‘■"L^ 



^L 



= p j- S ^cs + Cj. (2 ^ ^ ir. V^ ^ 



\y ‘ 

alt,o, W ein A (T' = sir p — 3 ^ 3^ , so that the entire equ'tion be- 
comte, 



W ^3^ 

- ^ — V rf = :^ 

ot L, 



^ In Ct 

CA S ir. Y ^ 



ein d" A r . 



^ P V s ^ 

The left side ie r.on-dIrr.ension&li.ied as before. Dividing the demon. nrtor 



V 



Y_ d y~ _ d_r _ y. 

tv df ° ' 



So finally there is the second tquaticn of mctionj 
A. 5 ^ - 1, (2 ! u - 3^^ - J, tin ^ r = C . 



The final equation of the longitudinal motion is obt'.ined bv equating the 
moments to the monert of momentum about the center of gravity cf the aircraft. 
..onsider.ng Figure 3 again, this ie. 



« g idfi—ill , M - 



3 y 






dt' 



W 2 2 

where, - i is the moment of inertia « nd i is the radius of jyration. Since 

3 y y 

M is a functiOTi of lift and drag which in turn are not affected by (T , neglect 
ing the influence of V, rnd uith the stick-fixed •’ttun tior, there is, 
o o . . . I 



W .2 d~( 0^ r ) d' ( 5N }r ) 

_ ^ V — 1 



6 7 



df 



jy 



dt 






= ^‘o ^ 



Sm . Sm . 



. ci jdi g” , 

oV oc^v d if 



This is seer to be the same as equ. tior '.3 i.ssumlng initi." 1 equilibrium. 



Putting the change in moment in coefficient form, 



V2 






V 5r-^ O L_ 5-^/ - 'bl t. — £2 — 

5v y 2 V N,d£ix 

^^dt^ 



P 



v£ 3/2 ^ 

2 ^ dt 



since, A c< V A r = = It * 

^n important point to note here is the Geman vee of the square root of 
the wing area as the additional length term when converting moments to coef- 
ficient fonn. The reasoning is that this quantity is more easily definable 
for the variegated types of wings now in existence whereas the man aerodynamic 
chord is somewhat nebulous in definition in the litorrture. Another reason 

1/2 

is the fact that the quantity r^S'^ ranges from C.8 to 1 1C and this facili- 
tates its approximation to unity where it appears. 



Considering the last term and non-dimensionali_ing the derivative bv lett- 

ing, dt V ’ 



Sc, 



M 






f 



vf. e.3/2 ^ 



dt 



= M pV 
dq ) 2 



o 3/2 ^ 



V: 



dt 



Then ncn-diraensionalizing the entire tern, 

p £ g3/2 -... y ..), _z_ i ^ - c p ^ 



y 



io< y ) 



Now ixon-dimencionalizing the left side of the equation, after dividing bv 



H 2 1_ d'PcK r ) ^ 

g V t; " 









= ((X 



r ) 



- 7: = • 



This equation nay be rewitten In the more concise form, 

^ ^ i Q i 1 

y :/ ^ .7 

There are then, the longitudinal equations of motion in nor.-dlmeneicnal fon 






f 



Ik 



A.k 



A. 5 



/S In D S In T 



“ - s <5-^ - f: s^> “ > « 



b In 



«'- =L ‘2 * s^rf'> “ • - “l ^ 



A. 6 ck - C^ ^ ^ " ^°io( ^ T2 ^* ■ S In V 

Q 1 ^1 0. i i 

7 7 7 7 



u = C 



These equations have been seen to compare term by tern, except for the 
}.ncLu3ion of variation of thrust with velocity, with these developed in a 
nore general sense. The assumptions and approximations have been cited. The 
dependent variables are u, o< and if'. 

A nore general case would be that of assuming, instead of the homogeneous 
set above, the equating of the quantities above to forcing functions of the 
elevator displacement, f . This would require, in addition, an equr.tion of 
motion involving moments about the elevator hinge line. This would be the 
case of stick-free stability eind motion. 

It is noted that a striking difference between these equations and those 
more fam-liar to most students is the form in which the stability derivatives 
have resulted. Since the problem is now ready for analysis of modes of motion, 
damping factors and associated desired results, the quest 'on arises as to why 
this form of the derivatives and how they are to be evaluated. J 

The usual procedvire is the separate evaluation of each derivative by 
similarity to previous determinations or by wind tijnnel or flight testing. It 
is seen that certain of the derivatives in the throe eqxxations above are in 
combination. For a special but usual case of the solution of the equations of 
motion further combinations of these derivatives appear. A discussion of both 
the general case and this special one of the solution foliows. 
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B, Solution of the Equations of Motion 

The most general approach to the solution of the equations of motion "be- 
comes quite involved in that the solution is carried to the determination of 
an equation for each variable composed of both the complimentary function and 
particular integral. See Reference 1. If the equations are restrictid to the 
stick-fixed case and the moticn is assumed resulting in the absence of an 
applied forcing function, the equations will be homogeneous and the solution 
resulting will be the complimentary fvnction only. 

The discussion here will be further restricted in that the equations in 
terns of the variables will not be the end result This would merely give 
the particular motion resulting from certain initial conditions. 

The characteristic equation which results from the solution of the dif- 
ferential equations being of the fom, x r x^e , will be analyzed for its 
roots and information relative to inspection of those roots. This will show 
the character of the motion as regards stabilltv, periodicity and demping. 

The natvire of the characteristic equation, or stability quartlc, will be 

Z 0 

set down from the determinant resulting from the solution, x = x^e Ex- 

perience has shown that the two sets of roots that form the solution are 
usually of such widely separate magnitude that they may be separated. The 
special case mentioned is baaed on this fact 

Since Scheubel's solutions employing the quasi-static stability concepts 
are based on the premise that the roots are separable, the general characteris 
tic equation will be set down without further comment and the special case 
will bo stressed in preparation for the solution by his quasi-static stability 
ccxicepts. The errors involved in this approximation to the roots, especially 
as it increases with angle of attack, is discussed et length in Psference 1 
It is shown that it is a good approximation for low angles of attack, i.e., 3° 
but becomes unacceptable for angles of attack near 10° especially for n«^utral 
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stability, 

L. The General Case. 

z'C 

Aestuning a solution of the form, x = x^e , vhere z ie a real or com- 
plea constant, the indicated operations on the variables are then carried out. 
The stability quartic rbsults from an expansion of the following determinant 
that follows from the fact that the equations are consistent If the detenn nant 
of their coefficients vanishes. 

Eeari’anging the equations by variables, 
u 



A .4 


. . ^ In D o S 

f In V D ^ 

< ° 


In 

In 


^Ju 




, « 


A. 5 


1 -=L <2 * 


d In C 

STFir) 


u 


• 


J' “C- sin If 

Ju, 


1 -c, • a - = c 


A. 6 

\ 


Sc 

in 






• 


* • O * 


^-c. -1 « -c^ ^ = c 

m . c. tn^ j 2 / 

q i w i 

7 7 


" h "in V 


^2 u 

7 




• 


^ O ^ 

rn 

’ *7 



The assumptions are made that sin s''- f is negligible since near horizontal 
flight is assumed, and that approximately, T = D . 

o , 2 

The determinant of the coefficients is then, using k“ = S/ i 

7 7 



I — 

„ /SlnD SlnTv 

i 

1 S In 

I ■ ^ l)''ln 

d In V ^y/" 

The qmrtic becomes, 




B.l 

vhere , 



4 

L. + 





a.k + a = 0 
1 o 
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„ ,1 p 1 ji D ^ In T x 

In V " h S V^ 

a< X ,2 . . „ In D ^ In T> 



-2 - ‘ ;;;r> S/>* S ‘rvf - - =- >'i> 



", y 



^ In C, 



* <'=L-S«) [°L<^*sTrr>] 



/„ ^ ^L(x V ,2 . , In D a In Ts 

^1 “ ^ " SI^Tv^ 



’C7( 



Sc 

ni 

a In V 



'J^ 

a In c. 



C (2 1. 

^ a In V 



-) 



a In Cl 

^2 a In V^ 

y/ ’ c 

^0(' 

a In C^ 



(C .c _£t)k2o..— 



^C/ 



“ c k" 

% - 



■5 a In C ' 

A (2 - nrv^) 



2 ^ ^ °L 

a = - cf (2 + ^ -- , - -■ ) 
o L 0 In V 



m 






he 

m 

a In V 






a In C^ 
<=I.(= ‘ njTv=> -‘ 



V/- 



The ueual taethods for the solution of quartics may ho employed to solve 
this equation as it stands. The constants have been arranged in the particu- 
lar fashion shovn for reasons to he mentioned later, 

2. The Case of Distinct Sets of Eoots, 

Experience has shovn that the usvial motion in flight consists of a long 
period, lov frequency, lightly damped motion, called the phygoid or flight 
pat*" oscillatlcn, and a short period, heavily damped motion, called the rotary 
oscillation. This suggests the factoring of the characteristic equation into 
two quadratic eqxiations. 

The length of the period of the rotary oscillation has been found to he 
at most of I -5 seconds while that of the phygoid usually is of the order of 
30,-60 eecands for most conventional planes. This suggests that, due to the 
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ah^t duration of tho rotary oscillation, the velocity ie essentially con- 
stant. Thus this node is governed b-' equations A. 5 and A. 6 with the terms in 
u elininatod. The solution for the phygoid can neither ignore the change in 
velocity or angle of attack. Eowever, consideration of equation A. 6 shows 
that, due to tho large oagnitude of the factor discussed in Appendix A, 
an approx.i nation may be made that the terms containing this item are large 
compared to the other terms. Solving this equation for ft' , the terms in (X 
in equatione P..h and A. 5 are then expressed by a tom in u. The phygoid ^alll 
be solved for, as regards the information mentioned previously, by considera- 
tion of the equations enclosed by thp dashed lines in the equation arro.y modi- 
fied by the approxinetion diecuosed above. The rotary oscillation will be con 
cemed by the equations enclosed by the dotted lines in the array, 
a. The Phugoid. 

Equation A. 6 becomes, nelgecting terms not involving yU , 



cr 




Equations A. 4 and A. 5 in coefficient form are now, assuming a solution' 



of the form, x = x e 

o 







8 In V 









In G, 



m 




5 In V 



0 









which gives the quadratic 



3.2 




I t., I 

o 



vhere 



1 



Sc 

m 

„ /d In D 0 In T^ „ In \ 

= S ^5"lTv ■ SHTv^ ■ ^~TT~^ 



n 



and, 



CK' 



a 



O d In 
, 2 /^ •’ L 



dc 

m 

a In V, 



■ * riKT-> * '=! '-= — > 

^ Cl ^ 



c^. 



C ?(2 



a In C_ 

L'i 

a In V ^ 



‘o< 



m 






ac 

m 

a In V 



C, (2 



c< 

a in 



a In V 



T^) 



The cuadi’atic has the roots, 



z = 



"" V 



^4-a 



“ -Pxl-ii, 

a. 



V/here the damping factor, is, = ■ “ the frequency or period factor is, 
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b. The Short Period OscLllrticn. 

iJet^loctlng the velocity terms, -and equation AA in its entirety, equations 
A. 5 and A. 6 hecome, in coefficient fOxTi, 

o -C. 



2 

' k ~ 

fliq J 



;c^Z C u 



= 0 



g 2 

v-hich gives, i- A^^Z + A^ =0 3-3 



« > ? 



where, A.^ = (j,. “ C. k ) an'’, A = (J C ) k t^.Cne root ie, z = 0, 

the others are, Z = <iam_ ing factor is, P^ = “ ~p; "the frequency 

factor, = Ja^ • f ^ 1 • 

'Before any further evaluation of the quantities determined up to this 



point, the quaei-static stability criteria are developed. 



20 



C. Quasi-Statlc Stability 

Before procoedirijj with the quaei-static stability development, it may be 
helpful If a reiteration be made of the Bpeciallzation assumed leading to 
this development compared with methods of a more ^^eneial natui'e as in references 
1 and 4. Both cf these references discuss the epecial cese of the apprcocinatlon 
to the roots resulting from the asBun^.ticn that the solution consists of two 
pairs of complex roots suffici&ntly separstod in magnitude to variant a factor- 
ing of the characterictics quartlc. This facile solution is usual in most 
cases but may not bo sound for unconver tional aircraft or even conventional , 
aircraft where additional facts must be considered. These are, flexing of the 
component parte such as wings and tail surfaces of the aircraft, which would 
entail additional degrees of freedom, large changes in the coefficients of the 
equations due to rapid depletion of a large fuel supplv and its effect on the 
derivatives. Many of these latter concideratione are undergoing exhaustive 
study at present and do not lend themselves to any concise gene ralix.at ions. 

From, say, a project engineer's point of view, especially in the initial stages 
of design, concern is giver, to information available from the most rapid and 
inexpensive data available. This Is especially true at present since specifi- 
cations, the military in pai-ticular, require employment of the finished design 
t>irovighout large ranges in altitude, speed and weight. The assumption of the 
presence of two quite widely separated modes, however, is sevnd for most conven- 
tlcrai aircraft and foi‘ unconventional aircraft at cei’tain phases of their 
flight histoiT-. The conditions that hold for the phugold were seen to be a 

Slav.' occlllaticn of long period and w.vc length and relatively light damping. 

« 

The rotary occilLetion is a short, heavily damped oscilla.tion. 

In ro^rd tc standard methode, the coefficients needed to determine the 

V 

constants for t'^e final solutioi.s cf the equations of motion oi'e obtained 
primarily from static vind tunnel tests. This is the rule for the values of 



equeticjnc A.l and A. 2 ae aeen in referencee 1, 3 The moment coeffici- 

ents are fo\aid by static etabilitv considerations for the variation of the 
EicJcier.t with (X . See reference 5* A dynamic \rind tunnel test is used to find 
C . See references 1, 4 and 6. 

The coefficients as they hare been developed in this paper will be evalu- 
ated in much the sane manner as regards those concerned with linear forces. 

✓ 

However, the method of evaluation of the logarithmic derivatives 'id. 11 bo il- 
lustrated. As for the moment derivatives, the detemination of C will not be 

m 

q 

discussed assuming that it may be found by the usual methods. The static long: 
tudinal stability criterion, , is the quanitity that the next few pages wi 

be concerned with. This quantity Dr. Scheubel determines in a quasi -static 
MB.J. Qauci-static in that it is developed from assumptions made on the d-mamic 
notion of the aircraft. The two concepts, quasi-static stability at equilibri- 
um and at ccnstaiit speed, follow from the fact that an aircraft has two dis- 
tinct phases of motion following a disturbance. 

Cuasi-static stability at equilibrium of forces is essentially a considera 
tion of the equilibrium conditions r'eached at a tine in the flight history 
following a disturbance at which the new velocity is obtained. This new veloc' 
ty is the initial value plus the incremental increase due to the perturbation. 
The quasi-static stability at constant speed development holds for the rels. - 
tively short period following the disturbance in which the velocity is assumed 
not to have reaches its final value. 

The d'.'namic response to an elevator deflection .a dicectod, so to speak, 
to glvo these two quacl-sts.tic concepts. The aircraft is initially in steady, 
nearly horizontal straight flight, starting from this steady stiTiight flight, 
we assume that the elevator angle, , bas been chango’d suddenly b a certa n 
amount, d £ , and we ask what will happen. 
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In tho first inctnnt, tho elevator displacement gives b chejage in longi- 
tudinal mcnent and a snail change in lift- This latter change ie so Ena.ll 
it is neglected. The change in nocient disturbing the equilibriun gives an 
angular velocity about the lateral axis, and, due to this, the angle of at- 
tack is changed too. A change in the angle of attack gives rise to a change 
of the lift coefficient and by it a change in lift. Soothe equilibrium of 
forces per jcndicular to the direction of flight ie dieturbed too. All these 
offecte hold during the initial portion of the period following tho distur- 
bance, i.e., during the first few eoconds. In addition, it is assumed that 
tho speed remains OEEcntially constant during this period. 

This then is the situation maintaining for quasi -static sta.bility at con- 
stant speed and considers the cheinge in force perpendiicular to the flight 
path caused by the initial curvature of the flight path and the change in mo- 
ments resulting from the dieburbancG. 

Eventually the change in elevator causing the change in angle of attack 
thereby the lift coefficient results in the aircraft attaining a new steady 
state, that is a new point of equilibrium of forces and moments at a certain 
speed, V dV, which is different from the initial one by dV. Cuasi-stctic 
stability at equilibrium conceme Itself with this portion of the flight his- 
tory. From experience it is known that it takes an aircraft an appreciable 
tine, normally several minutes, for the spaed to adjust itsoLf to a changed 
elevator displacement. 

1. Cuasi-etatic Stability at Equilibrium. 

Tho implication of an initial 
steady state of motion necessitates equi- 
librium of forces and of moments. From 
Figure 4 this is seen to be, ' 




'T 



Figure 4. 
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Ef: T-D-WsinJ^=0 

X 

Ef : L - W C 08 d*' = 0 

z 

Em : M = 0 

° ■ 



^Having etartud fi'on nearly horizontal flight the flight path will remain 
nearly horizontal for sufficiently snail changes so that the component of the 
weight perpendicular to the path remains alnost unchanged. From the second 
equation, that of equilibrium of forces perpendicular to the path of flight, 
then, 

d(L - W cos if*' ) = dL = 0, since, as before sin i“ = JT" = 0. 



now. 



dV 
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so that, f- 
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or, 
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g.dlnV.- 



then finally, 
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d CX 
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where the subscript E indicates an equilibrium consideration with the assump- 
tions implies above. This equation is neither a partial nor a total derivative 
from a mathematical point of view and caution must be exercised in handling it. 
If the new state of motion is a steady one, which it must be, 

dCj^ = M - = 0. The moment coefficient depends on 0^ , V, S , and the angu- 

d0 

lar velocity, ^ . Since, however, it is assumed that the flight path has no 
curvature at the nev/ steady state, end. Indeed, in actuality such would be the 
case, the angular velocity is zero. Thus the quantity, q, does not appear. 
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There ie then, 

“ 5^ ^ d In V = C 
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This then, is a nevr etahility requirement that applies to that phase of 
the longitudinal motion described and implied from the previous discussion. 
This requirement certainly is sound. For an aircraft which does not have thi 
quasi-static stability will be unstable in as much as for any disturbed state 
of notion, enforced by an elevator displacement, will show the tendency to 
move its elevator further in this direction, so increasing the deviation 
from its initial state. 

dC 

, The usual static stability criterion, is considered either, 

dCj^ dC^ dC^ dC^ 

- or, - (x - Xj, ). The quantity, , in either of the above ie 

L 

a constant throughout the usual range of normal flight, i.e, unstalled flight 
The value of x - 2 ^ can be found graphiclly from wind tunnel data of Cj^ vs. 

as explained in reference 5* This ie the more informative qmntity since 
a shift in the center of gravity is more readily determined throughout the 
flight history. 

The first tern in equation C.2 is very similar to the static stability 
criterion above. However for comparison, equation C.2 is written in its en- 
tirety as. 



I 

■ ~ — — I (x - X* 1 ) 

dps p dcK 1^ 'Z^ 



where Xj. |j, is a new neutral point that takes into account variation of the 
neutral point with velocity, i.e., compressibility effects. First, it is 
noted that variation of lift coefficient with angle of attack differs between 
the two. This is due to the fact that the ordinary determination of the change 
of lift coefficient with angle of attack is accomplished in wind tunnel tests 
where the stream velocity is kept constant. The results are then based on a 
constant dynamic pressure, whereas the present development is not. The dif- 
ference between the two is about l-3o/ o and is sumed negligible. The quantity 
Xjjlj. moves with respect to velocity, dynamic pressure, deflection of the fuse- 
lage, twist of the wings and horizontal stabilizer. A comparison with the 
movement of the usual static stability neutral point will be made only with 
respect to variation with velocity and dynamic pressure. Figvire 5 shows the 
variation with velocity of three different aspects of stability at sea level. 

I 

These are: (a) static stability only, (b) quasi-static stability assuming the 

zero lift moment coefficient, C„ , to be -0.02, and (c) quasi-static stabili- 
se 

ty with zero lift coefficient assumed zero. Figure 6 shows the comparison 

of the same quantities at altitude. Here the decrease in the stability margin 








H = lt-0,000 ft. 
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for qmsl-etatic stability is much lees due to the presence of the density in 

the dynamic pressure in the denominator of the additional term, i.e., in 
d In V 

E ■ 

2. Cuasi “Static Stability at Constant Speed. 

Kecalling the discussion of the conditions holding following the imposing 
of a small elevator deflection on the steady state, it was seen that in the firt 
small interval following the change the speed does not change appreciably. Thu« 
the quasi -static stability at equilibrium concept does not hold for this phase 
of the disturbed motion. It was further seen that the equilibrium of forces 
perpendicular to the direction of flight is disturbed too. The equation of mo- 
tion for this direction shows what happens. 

If a >forc 9 is exerted on a body moving in a certain direction, its path 
will be curved. See Figure If an 
aircraft, the ipoticn may become rather 

^ <A V- ^ 

complicated, for, due to the curvature ^ 



the component force of the weight per- 



pendicular to the path changes too. 



of the flight path in a vertical plane 




:i£_ I 



e 



Eestriction to small disturbances close 



\ 



to horizontal flight allows neglecting 



the changes in the weight component. 



Figure 7- 



For instance, a change in angle of flight path of 15° gives only 30/0 change. 
The centripetal acceleration is V . The equation of motion is therefore. 



and 



-V |t = d(L - W cos y ) = dL 

Q ^ V 










therefore, 
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and, multiplying by — - y - and noting that — = p — , there is 
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The term, — , may be neglected since is small compared to 
and, also, the usually large magnitude of^ makes it negligible compared to 
unity, j/bis is also the commcn approximation as seen by the neglecting of X 



and Z in references 1, 3 and 4. 



is, 



Thus there is the relation, , resulting from a considera- 

tion of the forces perpendicular to that of flight in the first instant follov 
ing a dlstui’bance. 

This curvature of the flight path or angular velocity has an influence on 
the equilibrium of moment about the lateral ajcis. So a consideration of the 
equilihrim of moments along this curved flight path involves the change in 
elevator deflection, the change in angle of attack and this angular velocity. 
The effect of velocity is omitted because of the assumptions of this phase of 



the motion. There is then, 

3c„ 3C., , 

np M ^ M ^ 



or, 
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This expressioi is similar to the one found for the equilibrium of a steady 
state of flight enforces by a certain elevator displn cement . It contains the 



stability tenri; 

C.h 



dcx 
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called "qua si -static stability at constant speed". A comparison of this quan- 
tity vlth the usual static stability is seen folloving a further evaluation. 

The angular velocity causes the angle of attack of the horizontal tail 
plain, iX to change by the s.mount, 



d « 
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dt ^h 



but, 
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therefore, 
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This change in tail angle of attack gives rise to an increase in lift coef- 
ficient of the horizontal tail plane, C.^ , and by it a decrease in moment about 

hi 

the lateral emcis, o 

53/2 

and, haviag notod its dependence on Angular velocity, there is, 
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dq 
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This quantity, , may bo determined by a dyoanic vind tunnel test de- 

q 

scribed in reference 1, or by use of the curved-flow teclinique referz’ed to in 
reference 6. 

dCj,^ do I 

The term, ^ j , added to the static stability is alunys negative, 

since Cj^ is a damping derivative. Thus this is an increase in static etabili- 

q 

ty and the neutral point corre spending to this quasi-static stability is be- 



. K 

I his apparent improvement over vind tunnel static stability ie seen to be 



hind t'^e one for static stability by about, — , where, K = C, 



larger for- an airex-aft of small ving loading at low altitude than for one of 
the high wing loading at high altitude. 



The additicnal tern added to the ueual ctatic et^bility neutral point 
looks suspiciously liko the ctick-fixed maneuver point, especially because of 
the statement of the situation and the eseumptione ro;ide. This is indeed the 
case and can be verified by conperiscn with the ncjieuver point found in ference 

k. 

Concerning the concepts of quasi -static st^ibilitv, it is seen that both 
are inprovenents over wind tunnel static stability. However, three facts should 
be borne in mind concerning their effect and their use. First, both contain 

i 

the usual stet'c stability plus an additional term Secondly, their application 
and use only hold for specific conditions perculiar to different phases of the 
flight histor:' Lastly, for a specific aircrc.ft they vary differently with like 
parameters, i.e., density, but all parametoi-s do not effect both 

Dr. Scheubel pointed that in Gemsany during the Inst war certain of their 
fighters showed a marked improvement in longitudinal stability over that pre- 
dicted by wind tiJinel tests at mediuai altutudes. Since these were of the k^G- 
500 mph class, this was explainable by the above quasl-stetic stability consid- 
erations. He said fuither, however, that pilots complained of poor stability 
at high spoece at high altutudes. In as much as both concepts still predict an 
improvement, greater for equilibrium and less for constant speed, the aeroneu- 
tical engineers v/ere greatly puzslod. Further study accounted for this by rea- 
sen of flexing of the wings, the tail pLane and the fuselage. The ac ountablli- 
ty of these phenomena in stability studies is very/ complex. For the longitudi- 
nal case a good approximation made be obtained from grouping their overall 
effect ac a shifting of the neutral point. Tho seriousness of their effect c<an 
be seen f.on mention of two simples aspects. First, upon the sudden npplicatior. 
of an elevator deflection or a change in speed as in a gust, the ©Lasticity of 
present designs at high speeds causes a twisting of the wing. This causes the 
angle of attack tu very spanvise resulting in a shifting of the center 'of 



3C 

r 

pressure. The same applies tc the horizontal tr.il surface. ' Secondly, £. 
shifting of neutral point vortically from flexing of the fuselage coupled vith 
the very largo drag forces at high speeds produces an appreclahle destabilizing 
nocicnt. 

That the problem is complex analytically is easily seen from merely a 
superficial study of the present day literature as, for instants, reference 7- 
Dr. Scheubel pointed out the experimental difficulties in a description of an 
attempt in Germany to determine the deflection of a horizontal stabilizer with , 
speed. This involved a complex optical system employing an arrangement of 
mirrors refloating a slit light source the movement of which was recorded by a 
motion picture camera. In addition to this the test aircraft was required to 
dive along a path determined by an ar.’ay of floodlights to insure perfect lin- 
earity of flight path.' 
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I). I valuation of the rerivativeE . 

Sone of the quantities to he evalu'jted are identic” 1 to those appearing 
in standard references. M-thods for the detemination of these will he men- 
tioned but not discussed in detail. 

Since all quantities appear in the general solution, the items \/lll he 
discussed in the order that they appear in the coefficients on page 11. 

. C_ . Static wind turnel tests. 



'(X 



M. 



lynamic vind tunnel test by the whirling arm or oeclll-ation me- 



thods mentioned in reference 1 or by the curved -flow technique referred to in 

t 

reference 6. 
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/ j.lct of ve. V can he obtained from static wind tunnel tests 
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Thus this varies f cm a value of 2 
at zaac M u,^ to about 3 fox M = 1.0. 

^ ki T ^ 

and ^ y ■ . Those aje cen- 
sidured to^.:ther hecauce they are found 
in nu^h tho same manner. From wind tun- 
nel tbcts ve. V are fouiia. Inline 
manufacturers usually submit dpt^ of T vs. V for specific installations. These 
derivatives are then determined as in Flexures 9 ?nd 1C. It is pointed cut that 
Fi^^e 10 is plotted for a constant of att^ick. 
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§~ ln ~ V ‘ ^ plo't of V8. V can be found from wind tunnel tests. The 



derivative is then found as shovn in F'guxe 11. 
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Figur'e 11. 
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E. Evaluaticm of the Longitudinal Modes Using Cuasi-Static Concepts. 

1. The Phugoid. « 

The damping factor, depends on the coefficient a., of equation B.2 and 
this quantity is, 



a. = ( 



B In D S in T 



D 'a In V ■ a In V 
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At relatively low velocity the approximate values of the derivatives are, 
= 2, = -1 (for a .let) end = 0. Thus the damping is. 



^1 " ""is 



At higher speeds the effect of the last quantity in will he felt. This 
will decrease the damping since it will overcome the increase due to the other 
quantities. 

The period factor, I^, depends on both the coefficients of equation’ 3.2. 
The coefficient, a^, using the quasi-static stability at equilibrium, equation 



C.2, becomes. 
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The period factor nov becomes. 
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Again, at low speeds the approximations can be made that. 
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So that. 
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The period is then, 



T = — 
P UJ 



and the wave length is, 
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2. The Short Period Osclllaticn. 



The damping factor is, = " 2 (^t " ^ )• 

^'p< w,, y 

Approximate values for these quantities are. 



Cl = 5.0, 






and, C k = -5 . 

q 



The damping factor is the, E^ = ■^' 5 - ("5) = "5 which indicates heavtr 

damping. 

The period or frequency factor was, 
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And, using the values given before, it is seen that the item in the parenthesis 
is zero. At least the value in the parenthesis is very small and it being to 
the second power, divided by four, in addition to the presence of the large 
factor,^^ , in the first term allows neglecting it in comparison with tv's firs 
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term. Thus there 



is nov, 
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From the foregoing it is seen that depends on: 

1 



l( Center of gravity location: 



yf 



X - 






2) Altitude: -~ . 

'^hl , c 

3( Wing loading: . 

4) Size: Hfp. 

5) Macs Concentration: 'y 

6) Velocity: ^ . 

The decrement is seen to be. 



.. <'=L^ ■ 
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Although this analysis is specialized and the modes of motion are of such 
frequency and damping that they do not present much concern in design, the ap- 
plication of the quasi-ctatic concepts has been shown. In regal’d to the im- 

•» 

portance of the modes, reference 8 points out that the importance of the phu- 
goid should not be m.nlmized too much. 
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The evaluaticn of tho coefficients of the non-fact ored qa'^rtic, eqa^tion 
B.l, page 11, in terras of the quasi-static stability concepts vlll not be car- 
ried out. The coefficients have been arranged in r form fer easy substitution 
of these quantities and a glance will show by vhet degree they can be simpli- 
fied by such a substitution. 

It is inp 02 -tant, however, that a criterion for the separation of the quar- 
tic into modes of slow and fast frequency be established. If the ratio of the 
imaginary pe.it of the roots descrying the phugoid to that of the short period 
oscillation is snail compared to unity, then the approximation is acceptable. 
From the simplified frequency factors used in the previous analysis, and, It 
is pointed out, the terms reV-ined are the most important ones, it can be seen 
t’"?t this criterion is primarily dependent upon the fhetors x - i., and^ in 
the nuraei’ator of This follcvs sinje the other terms, although the"' mv vary 

are relatively constant while the distance between the center of gra.vity and 
the neutral point may be anywhere from O.OO5 to O.C5. The large varlstion of 
the density factor is shoumi in Appendix A. 

The present-atlon thus far has dealt with the stick-fixed case only. The 
degree of freedom about the elevator hinge line is considered next. 






F. Introduction cf the ’ levntor Force. 

In Section P it ie recalled that the equr^ticn of loncitudinal motion due 
to the de^ec of freedoci about the elevator tos locked, i.e., it implied thst 
a change of elevator an^le irapreeeel and holi feet. Kov thie roetriction if 
removed and a rcfult obtained in the form of a inodifics tion to the original 
eolution of the ohert period oecillatior.. This is the stick-free case. It ie 
interesting to note here that the meet coimcn texts on this sublcct treat this 
caoe quite differently, reference *t- rules cut "nv modification of the phugoi'’ 
mode end solvco the remaining equations, dropping terms in u from an assumption 
of constant speed, for two short period modes. Cne of the Lotter is disregardec 
as the flapping cf the elevator about its hinge line, while the other is discus- 
sed extensively as the mode causing the "porpoising" notion, reference 1 
^ites off the freeing of elevators es a modification cf the modes by a reduc- 
tion of the tail effect by reason of the wing downvash and a possible contribu- 
tion to the velocity moment derivative, by reason of some constant moment 
due to elevator mass or a spring attachment or friction. 

The equation of motion resulting from this degree of freedom consists of 
equating the moments resulting from freeing the elevator to the inertial term, 
or moment of inertia of the elevator times the angulai’ ac eleiation. 



Thus, 
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Figure 12. 
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Fron the figvire It iB seen that a dovn (positive) elevatoi* deflection 
causes a nose doim (negative) aomsnt about the aircraft center of gravity. It 
was neglected for the stick-fixed case but will now be Included. This quantity 
a^ears in equation A. 6. Thera are effects in directions consistent with the 
other equations also but these are assumed small. 

It is also reasonable to assune that there is an increase in moment about 
the hinge line of the elevator due to the change in angle of attack of the tail. 
As is pointed out in reference this should’ be reduced to as near possible 
as possible. A typical fighter of the last World War had a value of C. 25 o/ c 
that of C, „ . C, may be made to be small by proeer design and many methods 
of reducing it are discussed in references k and 9- The incremental angle of 
attack of the tail itself is small for the portion of the flight histci^r con- 
sidered here, and even that is considerably reduced bv wing dovnwash depending 
on the horizontal tail's position. The C, term vlll be neglected here, but 
certainly should be included if design prevents its being made small /ny 
large value of C, will detract from the pilot's "feel" of the aircraft. 
Derivatives of the hinge moment with respect to rates of change of and S 

Xt 

arc nelgected also. 

As in reference h, any modlf Icr.tion of the phujeid from freeing the ele- 
vator will be resumed to bo of no consequence. Thus the improvicment to the 
sticV:-fixed approximation will only consider the short period osclllat on, 

vdiich, as mentioned before, assumes constant velocity. Thus this equation will 

✓ 

neglect all changes v:ith respect to velocity. ' 

The second term in Mg is the coriollls force which results from reference 
of the angular motion about the aircraft center of gravity to the center cf gra 
vltv of the elevator. The last term is the mass moment effect due to the angu- 
lar acclcration about the aircraft center of gravity, bo F.l becemes, letting 
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A = •- 5^1^ ~ and non-dimeneion? lizin^, 
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then F. 1 becooes, 

F.l t(l A)^-tB/^r^t(l^^)6< t^ S‘=0. 

Thic equation along vith equations A. 5 and A. 6 of Soctlcn A gives t^e svs- 

tem which shows t'^e modes of motion for the stick-free case. The complete S”S- 
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ten of equations, incl\iding A.k, should be considered for a rigorous £jial’'-8is 
of the modes of motion. However, as was diBcusGed pr-eviouely, the short peri- 
od oscillation presumes conetant velocity and this i£ ago,in aesvined. Thus thert 
is the following set of equations. 
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As Cuming a solution of the fozu, x = e , there is the determinrnt. 
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or, eliainating the root, z = 0, and letting, f,(z) = z + 2 (C, " “i_) 

i 1>C< 4 
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z) - C 



, Byuthere is. 



and, f2(z) = (l + la)(z + C 

F(z,t^) = f^(z) ^1^^^ ' “S/'^2^^0 “ 

If now the case is considered for the nass Aeneitj,^^^ , large, wliich is 

usually the case except for gliders, then the q\ae.ntity, t , essentially vaniehe. 

0 

Co, approximately then, f,(z) = C. The entire expression, F(z,t ), may be 
analyzed from the point of view of using the roots, z^, of f^(z) in the evalua- 
tion. 



The effect of freeing the elevators then will be the determination of the 
damping factor and frequency of the short period oscillation from a factor, ^z, 
to be determined. 

Beviewing the sitviation, the characteristic equ?.tion obtained from the de- 
teininant gives a fifth degree equation in z. Elimination of the null root 
gives F(z,t^) which is fourth degree in z. Beforence 4 then solves this equ'"- 
tion for two pairc of roots which give the two short peidod oscillr-tionE as 
pre V i ou s ly ment i oned . 
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There is then, in matheoetical terms, 

t = 0 giving, = 0 vhere. 
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It is seen that, since both ^ and C^ are negative, the real part, which 



indicateE damping, ie negative and thPe this motion is damped for this accep- 
tance of sign. However, the value of C might go positive. For instance, re- 

^8 

ference 9 shows that a horizontal tail design using an KACA 66-OO9 airfoil eecti 



with flap 0.30 chords in lengf-, an air aa.? of C.CO5 chords rnd with the hinge 



line at a positon of L.50 flsp chords, reduces the 



to iero between 



0< = -8 te o< = 5 , while C 






is positive. This reference also shows that 



^ nay go positive depending cai the use or non-use of ahovnbalar.ee. The 
use of a horn balejico, in general, appears to mr.ke C go slightly positive 
also, the magnitude being very saall 

Thus this mode of notion nvay be positively damped or possess moderate di- 
vergence. The fi'equency is given by. 
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where the ternn neglected is, (j^ - m ) , and is considered small oe in the 

previous diecusGion. The terms in the brncke-te are negative, the last term 
positive. The total va.lue of the quantity is not too large. 

It 'is difficult to strike a concrete conparlson of the mode of motion 
solved for here find thoce found in reference k. This solution does not solve 
the equation by a facto^i^ation but by the aesunr^tion of the quantity, t^, 
being large. In gcnora.l, the rcots solved for by Dr Scheubel are not f'e s-’me 
as these caueiIi^5 the porpoising motion described in reference 4 , since there it 
is indicated that the nature of the porpoising mode is a function of C. end 
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4h 



and C 



m 



It is not to be im- 



^ . Here, it is dependent upen ^ > ~ h ~ S 
plies that this is not such a type of motion, rather that it is not the same 
implicitly. 

Dr. Schcubel neglected the quantity 1 ^^^ , it is believed, duo to the fact 
that in Germany it was the practice to design the controls such that it was re- 
duced to a m nimun. He discussed theT h- jidling qu lities of " very large four- 

t 

engined flying boc.t In v/hich the stick forces per g wore very cmrll without 
po\7cr boost, tbc exact values of ^/hick are not remembered. In regard to 
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lateral control, hovcvor, it is remenberod that he stated that one person 

covld EKilntain straight and level flight with both engines dead on one side 

with one foot on the rurllcr control, again vithout power boosted controls. 

Inclucion of the 0, torn would plrce it in the solution of f,(z) so 

‘■‘cK ^ 

that it would detirxt or add to both the damping and frequency terns as it was 

negative or oositive in the rrtio <;f Furthermore, in os much as down- 

wash would multiply the value of the variable, 0 ( , which rffects it by a fee- 
tor less than unltv, this would fui'ther decrease the ratio itself. 
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G. Effect of Elcv«.t Inpulse Typo Forcing Function. 

Ecitei'citing, it is recalled thst the stick-fixed c; se involved the impor- 
tance of en clov- tor dlspleceaor.t or a guat one the ensuing motion v- e nnnl”-- 
zed in tiTC light cf the elevrt,r being locked in position. The etick-fr-ee case 
allowed the elev't>r to float freely abo\:t its hinge line. Furthermore, the 
equations of notion were ar.runed honogeneoue. /'n additional c.'^se Is now con- 
sidered in which the equations involve a forcing function ind further this 
forcing function is on el,vf tor impulse, an elevator def locticn of short dura- 
tion. This case is analogous to the incorporation in the system of ■'n on-off 
type autopilot. 



An elevator impulse as shown in Figure 13 le imposed, 
given the cyrahcl, I ^ , and is. 

If = / r d>r 
0 

or. 



ir • = /J s dt 



This impulse is 



^ The first phase folloving the impulse, 
that of the short-period oscillation, is 
cor.sidered and thus there are the equations, 

A. 5 r ■ sir . 



Figure 13- 

= ST 
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^.6 - m^<k - ciqyiACX = yU S { 'o ) . 

The effect of the in,.ul6e on the eqxr.tion of motion perpend icuLar tc the 
wind axis is negligihle for conventional nircrrr.ft since the only force created 
is a lift force, -i/hich, due to the short duration of the Impulse, does not not 
Icng enough to be affective. This is not so fo tnilleeo r Ircr ft, however, 
since for ruch aircraft the deflection of the elov'^tor, w'hich is '•ctually a. 
portion of the main lifing foil trailing edge, changes the wing angle of '•tVek 



and causer nn ap,;reciable effect in this direction. 



he 



The effect of the impulse on the equation of monente causes rn inltinl 
angular velocity which vill he designated , hut no anguL”'r displacement. 

If there is no disturbance in flight path angle, which ie nearly correct for 
conventional aircraft, this initial angular velocity will consist essentially 
of solely a time rate of change of angle of attack. This la go becr.use of the 
mass distribution of a conventional aircraft compared to, s®y, a talllesB air- 
craft. Thus Initially, at = 0; = 0; c>(^ =0; o( = u>o 

Solving A. 5 for , eubstituting in A. 6, there is then for the short per- 
iod oscillation. 
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* (B= . Jv'l = - 



The solution for 0< le. 



0< =» e (A^ cos A sin )• 



Imposing the first initial condition gives, A^^ = 0. Differentiating then, 

. ^ E.. ^ 

o( = = e ^ ^^1 *1 ^ 1 ^ ' ®i ^ 1 ^^ • 
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The second initial condition gives, B. = . The Principle of Impulse and 

•^1 

Momentum states that the change in angular momentum about any axis equals the 



angular imt/ulses of the applied forces about the axis, or, 
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:: cv^^luatoc the constant, 3^, so th^ t consequently, 



From equation A. 5 there ie, ^ = c 
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Since at, 'Z = 0 , JT == C there is flivally, 
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This initial phase of the response to an impulse type forcing fvmction is 
depicted qualitatively in Figure ll»-. In equation G.l the vr"'luee of and 
ere the same as those found in Section 
F., page 3^- Thus the transient is 
heavily dam:v2d and typical values show 
it to have reached tho steady state in 
one -ha If periois. 

Examining further the steady state 
of G. 1, 
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If tho de- 
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nominator ic email, the roc ponce is very sensitive to the impulse. Dr Echeu- 

bol cited an example of neutral ctahility where the steady state then 'bocomoB, 
C I 5- 
X- = ^ 

^ 1 eince, x = x^j 
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Lsing values of. 



V — hyO mph^ o = 27c ft. i I5 = -1 ; t_^ = 0.2 sec. ; u = 
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the response is, 
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If • t = -c. / X 1 X ^ X - ^?.C035 = 5-7° 
5 a lt> . 5 



which is appreciable fo.'’ such n smll iiapulee. It ia pointed out that thie 

=^P(h 

value of elevator cffcctivonese, — , corroepon. s to rn elevator area, to to- 

tal horizontal tail area ratio of C.6 according to an empiricallv determined 
curve in reference h. 

The eccond phase of thv motion following the impulse ie that of the flight 
path oscillation. Thus equations A. 4 and A. 5 are used. These should be modi- 
fied ae before for the variable, c\ . Towevor, it will be o-ssumed that the de- 
rivativo, ^ , is approximately zero and will be neglected. 
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Fran the pact analysic of the two modes, the rel- tlve magnitudes of t'leir 
periods, and the magnitudes of their damping factors, it can be assumed, that 
the response t^, the impulse of the rotary motion liac reached the stdady state 
condition by the tine the flight pa.th response takes effect, to the Initial 
conditions are, at, % - = y* , and u = 0, Ck.n 'f' •=. o , 

Again the fcrcin^ function creates no nccountrble effect cither along the 
wind axis or perpendicular to it. Thus a ,,ood approxim. tioii is to verk with 



honogcncous cassations, do. 
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A solution for t ie, ^ = 
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The inclination of the flight ,>f th, , nr.y be f'cfine^l, using Fljuro 1^, 
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This rceult trsT-y be obtoined by a ciitjhtb' different appra-ch.' ^oneiderinc equa- 
tl' i! /■ , 5 obove, 
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Differentiating G.2, 
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and, neglecting energy loss due to draping, 
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kco that, l(h 1 — ) = G or, dh=- — — =-“U hnd, finall- 
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which ie s.e before. 



Using the seme figures ac before, the mrrxicrum gain in altiti^de ver- ap- 
proximate l,v b^coracc, 

_h = ;< 0.0995 = 9 % feet 
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Diecucsicn and Jonclucions 

The quasi-Etatic etabilltv ccjncct^ts baeed on equilibrium of forcee Pnd no- 
raents at varicuc phaccc of the fli^^ht hictor^'^ of ci\ aircraft follc^.in£' a dis- 
turbance appear to bo sound for an accepb-'nae of the description of the flight 
hietoi’y stated earlier, rioglectint^ adc’itionrl effects such "'e elastic It", these 
concepts account for the ap'parent increase in ct^hilit" due to compreeeibillt" 
effects associated vith the aircraft drga, lift an' moment If t^ie v lues of 
these quasi-static stability criteria could be "scert ined individually and in 
their entiret", ,^reEu.mably they vculd include intei-forence effecte. As it is, 
no American litei-atuiT known to the vn-iter discusses lonjituiin'’ 1 st'’bilit" 
comparisene at lov; and hijh speeds Dr. Scheubel rl^scusets the comparison wish 
the results sho\’?r in F djuros 5 f'-nd 6. This vve apparentl" b'^sed cn results of 
flight tests or- tho inclusion of the additional factors in the qw? ei-etati’c 
criteria anr.lytically dctermln d. This 1 tter would neceesaril" neglect in- 
terference effects. ' 

For the special c<ace of separation of the motion into the phugoid rjid ehor 
period modes good agi'ceaent v;ith standard developments vase had neglecting t'^e 
additional factors mentioned above. 

The modification of the phugoii modo due to f-xeing t^-e clev'tore \7re ne- 
glected as is usually done. The resulting equation vjoe anoly. ed for the epe- 
cial case of the aircraft maao dencity factor of 1 rgc magnitude. This assum- 
ption vas used to approxlcatc only one po.ir of the complex con'ugote roots 
Comparison with standard references shoved this motion to he approiimatel" the 
"porpo/Sing" node although not implxcitlv. 

Cf some intercet ic tho Gem.'’n use of the quantity 3“^, square roct of the 
wing area, ac the non-iinensionalizlng length term in conversion of forces and 
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laonontE to cocff icicntr. Actvc.ll’'^, Dr. Schou'oel did not point this out ce 
having vide Gerrm ueage, but rfthir ne being used et raxTri8t'’dt Techn^C"'! Iij- 
otitute. It m.ght be pointed out thr t the concept of rspect rrtio would be 
loct \>y uce of t''l£ Icng'-h. This night be bettor than the ncbuloue definition 
of noan aerodynamic cl;orrl present in current American aerod-mamJe literature 
For inctance, reference 1+, suppooedly a fairly rccredited text, defines moan 
eerodynanic chord ae 1) an average cherd^ for a finite wing of varying section, 
which, when nvltiplied by the aveinge section moment coefficient/ dynrnic pres- 
sure and wing area, gives the moment for the entire wingj 2) the chord of an 
imaginary airfoil which throughout the normal flight iTinge has the eemo force 
vectors as the three-dimenrional wing and can be determined from vind. tunnel 
teste or a gi’aohical analysis which asaunse a special shape to the wing plan, 
and 3) & chord on which all th(? forces and moments acting on the surface c"n 

be i-epresented as acting and on which the pitching moment coefficient is in- 
variant vrith lift coefficient. Dr. Scheubel pointed out that the special rea- 
son foi’ the use of this quiantity was the fact that the ratio of the square root 
of the wing area to thu tail length was found to range befv^cn C.C and 1 25 and 
the average value found for seme ninety airciaft irae 1.01 The advantage of 
ap;.roximating this ratio to xmity resulted in the adoption of this length. 

The extension of such a method to the lateral or asvranetric case is d.is- 
cuseod in gcnci'al tenne in Appendix B. 

The development for the longitudinal case as c. whole does not consider all 
the factors bearing on the problem by any means. It was pointed out by Dr 
Scheubel that this was a good apprcxlmntion from the prolect engineer's view- 
point and must necessarily bo improved as the design develops. For instance, 
tho effects of dcwi^/ash and the influence of the povei’plnnt have a Large ef- 
fect on the damping term. In any case it is a slightlv different view of the 



problem and its voirhrcsts on that. 
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Arrn-Tix ^ 



The nacs nunhci-j^yu , as mentioned in reference 1 'wr.e first exic;eated hy 
Gla.uert for ucc in the non-iimcnsionalizinc of the equations of motion That 
it also follcvt from dimensionc.! analysis, e.^-, the Pi Theorem, is well known 
Due to its extensive use in the litei'cture, a discussion of It here is onl''- vo.r* 
rantod hy the slightly diffcx-ent definition of it in this paper and the fact 
that considerations of its raasnltude fox' uoe in the approximations in thepape^r 
necessitate rccallinj vrhat factors bear on its magnitude. 

The diffei'c-ncc in definition lies in the use of the mean wing chord as 
the repreoentrtive length in roferonces 1 and 4. Furthermore, the masE or den- 
sity number ucod in this paper is tv/ice that used in these references 
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from which it is seen to depend on altitude, ^ , on iring loading, V//S, and 
on weight, l/W. The tendencies with these parameters varying are as follows; 

a. Altitude - As altitude increases, increases. An altitude increase from 
see. level to 40,000 feet causes an Increase in of from 1 to 4 

b. Win^' Load'ing - Since it is to the 3/2 powoi', an increase of wing loading 
of 1:2 gives on increase in^^ of 1 : 2 . 83 . 

i 

c. Weight - Ax, weight Incroases,^^^ decreases as (l/W)^. In this case a 1:2 
increase in weight give a 1:0-7 decrease in yA . Tj’-plcal values of^ ore: 

Lergc gliders at low altitude : = 5 - 10 

Snail fighters at high altitude: y/. = 5 OC 
Large ccrxicrcial pL'’.nes; yui = 50'1C>0 
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APPlliDIX B 



Biscucsir.n of the Extension of (.ucsi-LV tic St. hility to the Lateral 



or Asy.imet?lc Motion 



The possibility of the oxtensicn of the quasi-static stability concepts 
to t’-^c aajr;pxictric motion of an aircraft seems to be remote insofar ra the 
sane line of reasoning be followed. It will be recalled that in both quaei- 
staiic stability/- determinations consideration made of forces in a single 
linear direction and of momenta about a point in one plane. In effect the 
ratio of changes in parameters affecting both sj/etenia lead to the quasi -sta- 
tic stability concepts. Cf fundamental iaport-ance wric the fact that consider- 
ation of those changes did net affect the additional linear force equation 
Cuch a correspondence does not hold in the lateral case. . 

Aesuming a conventional aircraft it has been found that the modes of mo- 
tion in the lateral case are four in number. Following the same line of rea- 
soning as was done in the longitudinal case, a disturbance could, all factors 

V 

ccncidored, cause a linear displacement laterally, i.e., along the y-axis, a 
rolling displacement, and a yawing displacement. References h, 10 and 11 show 
that the solution of the lateral equations results in a characteristic quar- 
tic which gives roots X ^ 1+ ’ solving for these roots 

all neglect the presence of products of inertia, that is, they assume that the 
principal longitudinal axis of the' aircraft in is line with the flight path, 
rnd assume level flight or tan 0 = 0. Tl:e moment composed of three components 
given by: 
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ty cf rcll, rn'^ n< 3 ul?r velocity of yru roepectlvely. Ihe coefficients 

ojid rre F o^endent en initlc.1 conditions, on’v four '.icing ^ r'^ • tr- I'y , the 
reminder are detcinlnod fion these. These individuc.1 motions urill he con- 
sidered sCficretcly. 

Considerin’ f'c first group, the root, "is err 11 rnd ueir lly positive 

giving light divergence. This is the so-called spiral motion end either ?. scrl 
dieturhanco vrtiich io j,uio sideslip, roll, or y; v will develop cn interdepen- 
dence of hoth rolling end yawing moments on heth sideslip angle end yeving 
velocity. Following q,uaei-atatic reasoning, since the period is long, t’-ls 
metier indicates t’ e ertablishment of a unique rcL' tionehip "between initial 
disturbance increments and resulting forces and/or mucierts in equlir^rium 
after a long period of timo. This does not follow since either ailoron, rud- 
der cr sideslip, due to a gust, sav, my cause the motion and all freots of 
lateral motion aie involved in the resulting moticn. 

The second group ic characterized by a large root, A ^ ^usually negrtivo. 

Bg ard are smll in compel ieon ^.ith so f-at the motion is heavil;’ damped 
and is escentinlly pure roll. To known reference cor.eideic this motion, of on- 
importance. Fefcrcnce 12 w'hich includes the effect cf finite product of In- 
ci-tia and initial inclination of the flight p.^th cxcliJes "ny mention of It leo. 
Pefcrcncc 13, which ic an investigation of unsrticfact'.iy 1 tor'll at- b lilt- 
on an actual cur ent cor.br.t aircraft, does net mention it. This evidence 
ccupled with t'^c fact that the motion involves only one degree of freedom 
ceMcalnly door, net lend itself to any quasi -e-b; tic stability criterion. 

The tbirO, group ic concemod with the cll-lnport-ant oscillator'^ motion 
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For fixed controls, this ccnciote of both a rolling action and a yawing notion. 
It tcc cm be induced by any one or any conbination of the pocsible lateral 
dictui’bancec. 

Fron the foregoing it io concluded thntunique representations in the quas 
static acnce cannot be deduced for the asyoLietric motion. It appears from the 
conclucions of reference 13 that a poccibilit'’’ n ght be a consideration of the 
degree of fredom abov.t the rudder hinge line. In the investigation outlined 
in t’ is report, r^:^dder“f ixed considerations were ruled cut in that actual 
flight tccts showed good agreonent with predicted oscillations up to Mach num- 
bers corresponding to a speed of a?h (density altitude at this speed •vra.s 
not specified). Insofar as a gradui'l deterioration in damping with speed 
result^^ng in ruddcr-freo tccts, a variation of the rudder hinge moments, 
both C, and C, , with Mach number is Inferred. It is also stated in 
this report that this trend has been observed in high Mach number teste of 
ether contraol surfaces. Amplification on these latter raentionod tests is 
not riade. IhfO facte must be pointed out here, ^ however . First, the above men- 
tioned possibility is not in the quasi-static sense of Dr. Scheubel's approach 
Secondly, reference I 3 based its computed values at high speed on an extra- 
polation of a determination at an airspeed of 210 miles per hpur assuming a 
constant value of the stability derivatives throughout the speed range 

A fx^rtihor coraplication in the establishment of quasi -static criteria fer 
the lateral case are the inclusion or exclusion of the- effects of products 
cf inertia. These effects are disciujsod in detail in reference ll4- and men- 
tioned in reference 12 for a specific aircraft conf i(_ur''tion. The effects 
of f^eir inclusion was considered negligible in reference I 3 , but this con- 
clusion must be made weighed in the li^'ht of the remarks made above. 
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